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We investigate the consequences of non-minimal gravitational coupling to matter and study how 
it differs from the case of minimal coupling by choosing certain simple forms for the nature of 
coupling, The values of the parameters are specified at z = (present epoch) and the equations are 
evolved backwards to calculate the evolution of cosmological parameters. We find that the Hubble 
parameter evolves more slowly in non-minimal coupling case as compared to the minimal coupling 
case. In both the cases, the universe accelerates around present time, and enters the decelerating 
regime in the past. Using the latest Union2 dataset for supernova Type la observations as well as the 
data for baryon acoustic oscillation (BAO) from SDSS observations, we constraint the parameters 
of Linder exponential model in the two different approaches. We find that there is a upper bound on 
model parameter in minimal coupling. But for non-minimal coupling case, there is range of allowed 
values for the model parameter. 
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I. INTRODUCTION 



The gravity action hence becomes. 



It has now become fairly well established that the Uni- 
verse is undergoing an accelerated expansion in recent 
times (z < 4). Observational evidence for this mainly 
comes from supernovae la m and Cosmic Microwave 
Background Anisotropics Q- Large Scale Structure for- 
mation %?] , Baryon Oscillations [4l| and Weak Lensing [5] 
also suggest such an accelerated expansion of the Uni- 
verse. One of the most challenging problems of modern 
cosmology is to identify the cause of this late time acceler- 
ation. Many theoretical approaches have been employed 
to explain the phenomenon of late time cosmic accelera- 
tion. A positive cosmological constant can lead to accel- 
erated expansion of the universe but it is plagued by the 
fine tuning problem 6]. The cosmological constant may 
either be interpreted geometrically as modifying the left 
hand side of Einstein's equation or as a kinematic term on 
the right hand side with the equation of state parameter 
w = — 1. The second approach can further be generalized 
by considering a source term with an equation of state 
parameter w < —1/3. Such kinds of source terms have 
collectively come to be known as Dark Energy. Various 
scalar field models of dark energy have been considered 
in literature [7h23] . As an alternative to dark energy as a 
source for the accelerated expansion, modification of the 
gravity part of the action has also been attempted [2J] . 
In these models, in addition to the scalar curvature R in 
the gravity lagrangian there is an additional term f{R). 
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S=^ j d^x^—g[R + f{R)], 



(1) 



However, in such models matter and gravity are still min- 
imally coupled. Despite the significant literature on such 
f{R) models 2J], another interesting possibility which 



has not received due attention until recent times is a non 
-minimum coupling between the scalar curvature and the 
matter lagrangian density [23 • 

In this paper we study the evolution of Hubble parame- 
ter in minimal and non-minimal coupling between scalar 
curvature and the matter lagrangian density. We use a 
form for /(i?) model which assumes the following expo- 
nential form proposed by Linder [26|: 



f{R) = -C[l- exp(-i?/i?o)] 



(2) 



with C being the model parameter and Rq is the present 
day curvature scale. We also attempt to place obser- 
vational constraints on the parameters of this model in 
both minimal and non-minimal coupling of scalar curva- 
ture with matter lagrangian density. We find that there 
is an upper bound on model parameter C in the min- 
imal coupling case. For non-minimal coupling between 
matter and gravity, there is a range of values which the 
parameter C is allowed to take. These bounds depend 
on the present day value of q^ that we use. The paper 
is organised as follows: In section [H] we introduce the 
most general action for modified gravity. The equations 
of motions corresponding to this action are solved nu- 
merically for both minimal and non-minimal coupling of 
scalar curvature with matter. We investigate the obser- 
vational constraints on the model parameter in section 
mil In section HVl we summarize the results. 



II. F{R) GRAVITY MODELS 



We start with the general action for modified gravity 
where the curvature is in general coupled with the matter 
lagrangian: 



S = d xy/—g 



h{R) + [1 + f2{R)]C,, 



(3) 



where fi{R) and f2{R) are the arbitrary functions of 
Ricci scalar R and Cm is the lagrangian density for mat- 
ter which we will assume to be non-relativistic. We as- 
sume the natural unit with the speed of light taken to be 
unity in our calculations. The standard Einstein-Hilbert 
action is recovered with /2 (i?) — and /i {R) = -^ where 
K^ — SnG. The standard f{R) gravity class of models 
are recovered with f2{R) — and /i(i?) — R + f{R), 
where f{R) is an arbitrary function of R. In the latter 
case the pure gravity action has a non-minimal coupling 
while the matter is still minimally coupled to gravity. 
In what follows, we shall consider the cosmological evo- 
lutions and their observational constraints for modified 
gravity models in both cases, namely, one in which the 
curvature is coupled minimally as well as the one in which 
the curvature is non-minimally coupled with the matter 
lagrangian density. 



A. Minimally Coupled f(R) gravity models 

As mentioned earlier, we recover the standard f{R) 
gravity (i.e. the one with a minimal coupling of matter 
with gravity) if we have f2{R) = and fi=R+ f{R) 
in the action defined in equation ([3|). Now varying the 
action given in equation (|3)) with respect to the metric 
tensor g^^, we get the modified Einstein equation: 



f 
G^u + iRRf^u - (^ - n/flJSp^ - V,,VufR = i^^Tf,^: (4) 

where fn = -^ and fjm = ^. Assuming a flat 

Friedmann-Robertson- Walker spacetime with a scale fac- 
tor a{t): 



today. It will be convenient if we express the cosmological 
quantities involved in dimensionless units. Hence, we 
define the following dimensionless quantities: 
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(10) 



fxx = -T^ — R0-TB2 ~ RoIrR (11) 
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a= 772 = 6(1 -go), 



(12) 



where Rq is the present curvature scalar and go is the 

present day deceleration parameter. The evolution of R 
with redshift, z is given by the equation. 



dR , dH , , , dli , n . . d H . , . 

-=6(3i7--(l + .)(-)^-(l + .)i7^). (13) 

With these definitions, one can write equation © in 
terms of the redshift z, as. 



dh 



K' + il- qo)f + /,((1 + z)h— - h^) + 

dz 



a + ^y 

1 - go 



1 ,dh^^ d^h 3 dh 

dz"^ I + z dz 



h'UAT(^r + ^- 



'' h dz' 



) = rtmii + zf. 

(14) 



Instead of working with a general f{R), it will be more 
useful if we use a specific form for f{R). A simple form is 
case of "Exponential Gravity" proposed by Linder [26| . 
In this case the form of f{R) is given by, 



/(i?) = -C[l-exp(-i?/i?o)] 



(15) 



This can be expressed in the dimensionless form by defin- 
ing a dimensionless constant c = C/Rq. We then have 



ds^ = -dt^ + a{t){dx^ + dy^ + dz^), (5) 

the 0-0 component of the modified Einstein equation ^ 
becomes 

H' + ^-fRiHH' + H') + H'fnj,R' = ^^, (6) 
D a-^ 

where / is with respect to In a, flm is the matter energy 

density parameter today and Hq is the Hubble parameter 



/(x) = -c [1 - exp(-a;)] , 



(16) 



With this choice of f{x), we now solve equation P^ to 
find H{z). There are a number of investigations where 
this equation has been solved and the model is sub- 
sequently constrained by observational data. In all of 
these works, it is assumed that the universe behaves as 
a ACDM model in the past and subsequently deviates 
from that behaviour [28| . In this way, one sets the initial 
conditions for h{z) and dh(z)/dz, assuming the model is 
close to ACDM in the past. One problem one usually 
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FIG. 1: Behaviour of Hubble parameter as a function of redshift. ilm = 0.25 and qo = —0.55 for both figures. In the left panel, 
c — 1.1, 1.2, 1.4 from left to right whereas in the right panel c — 1.6, 3, 10, 50 from right to left. 



faces in this approach, is the epoch at which to set the 
initial conditions. Depending upon the redshift at which 
one fixes the initial conditions, one may or may not get 
well-behaved solutions without any pathology. In other 
words, one has to fine tune the initial conditions so as to 
get regular solutions. Although this always happens in 
most of the relevant modified gravity models discussed 
in the literature, it has not been sufRciently emphasized 
to the best of our knowledge. To circumvent this prob- 
lem, we take a different approach to solve the equation 
P^ . We set the initial conditions at present epoch, i.e 
at z = 0. In equation (fT4|) . h{z = 0) = 1 identically. 
Now for the second initial condition, dh/dz \z=o^ one can 
write dh/dz{z = 0) = 1 + go, where go is the present day 
deceleration parameter, as mentioned before. Hence, the 
second initial condition is directly dependent on the de- 
celaration parameter at present and this will be one of 
the parameters in our model. Hence, we have three pa- 
rameters in our model, i.e., c, which is the parameter 
in the Exponential Gravity model, and two cosmological 
parameters, Q.m and go- Using these conditions for to- 
day's epoch, we evolve our system from present day (i.e., 
at z = 0) backwards in time (i.e., for increasing z). We 
want to stress that, in this approach, there is no extra as- 
sumption, or fine tuning of the initial conditions in order 
to solve the evolution equation (14). One of the initial 
conditions, h{z = 0) is fixed to 1 by definition and the 
other one is related to goi which we shall constraint by 
observational data. 

We aim to investigate if the cosmological evolution is 
well behaved as we go to earlier times, or whether they 
are plagued by singularities at high redshifts. We further 
aim to study the role of the parameter c in this context. 
As we discuss below, one indeed gets regular solutions 
upto any higher redshift for a certain range of values for 
the parameter c. For the rest of two parameters, we vary 
Q.m between 0.25 and .0.35 and the present deceleration 
parameter go between —0.9 and —0.55. 

We first investigate the behaviour of the normalized 
Hubble parameter {h{z) = H{z)/Hq) as function of red- 
shift for c < 1.6 as shown in figure 1. We have used 



three values for c, 1.1, 1.2 and 1.4 and we have assumed 
fim — 0.25 and go = —0.55. It clearly shows that sin- 
gularity occurs at different redshifts for different values 
of c. Even if we vary U,m and go in the range mentioned 
above, the overall behaviour does not change in the sense 
that there is always a singularity at some redshift. How- 
ever, the behaviour significantly changes once we assume 
c > 1.6, as we show in figure 2. In figure 2 we have 
again plotted h{z) as a function of z but with the value 
of c greater than 1.6. We evolve the system for redshift as 
high as z = 1000, and the behaviour of h{z) remains well- 
behaved. As we mentioned earlier, although these plots 
are for specific choices for ^Im and go, the behaviour of 
h{z) is similar for any value of these two parameters in 
the range mentioned above, i.e., 0.25 < Q,m < 0.35 and 
—0.9 < go < —0.55. So we can conclude that for the 
model parameter c > 1.6, the model is regular upto any 
higher redshift. 

We also study the behaviour of the deceleration pa- 
rameter g(z) as a function of redshift. In terms of the 
redshift, go is given by. 



g(z) = (l + z) 



H{zO) 



- 1 



(17) 



The behaviour of g(z) is shown in figure 2 for different 
values of the model parameter c assuming fl„i = 0.25 
and go = —0.55. It shows that in all cases, the universe 
has an accelerated phase at present and as we go back 
it smoothly enters the decelerating phase. The lower the 
value of c, the universe enters the accelerated phase ear- 
lier. 



B. Models with Non-Minimal matter-curvature 
coupling 

Here we extend the f{R) gravity models assuming a 
non-minimal coupling between the matter lagrangian Cm 
and the scalar curvature. We assume the action to be. 





FIG. 2: Behaviour of the deceleration parameter q{z) as a 
function of redshift. The hues are for c — 2,5, 10 from bottom 
to top. we assume ilm = 0.25 and go ~ —0.55. 
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FIG. 3: Behaviour of the Hubble parameter q{z) as a function 
of redshift. The lines are for c — 1.5,6,15,50 from top to 
bottom, we assume Qrn =0.25 and qo = —0.55. 
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,^R+il + f{R))Cr, 



(18) 



Minimizing the above action with respect to g^^ gives 
the modified version of the Einstein's Equation: 



FIG. 4: Behaviour of the Hubble parameter q{z) as a function 
of redshift for minimally coupled (solid line) and the non- 
minimally coupled case (dashed line). We assume Qm = 0.25 
and go = —0.55, c = 5. 



which results exchange of energy between the matter and 
the scalar degrees of freedom present due to the f{R) 
gravity model. This exchange of energy is one important 
feature of this non-minimally couple f{R) gravity model. 
But once we assume the prefect fluid form for our mat- 
ter energy momentum tensor (which is consistent with a 
homogeneous and isotropic universe) together with the 
form for the lagrangian density £,„ = —Pm ifim being 
the matter energy density) [27| , one can explicitly show 
that putting v — Q in the above equation, i.e for the 
energy density conservation equation, one gets the usual 
equation, p™ + 3Hp,n = 0. With this together with the 
metric given by ([S]), one can write the 0-0 component of 
the equation ([TO)) : 



H<l>+-fi~ 



(21) 
where dot represents differentiation with respect to time. 

Further, comparing equation [3] and I18[ /i is identified 
with _R/k^ and /2 with /(i?) given in equation[T5l Chang- 
ing the variable to redshift z, and expressing everything 
in terms of the dimensional variables defined in ([7]) , one 
can now get from equation (pij) as 



0i?M- - i^Rgt^u = k'(1 + f{R))T,,u + (V^V, - .g,,.n)0, 

(19) 
where (j) — l + 2K^Cm.-j^- The Bianchi identity G^u;u = 



gives 



VpT^. = Y^tII) (.9^-^™ - T^.)V^i?, (20) 

which implies the non-conservation of the matter energy 
momentum tensor. This is due to the non-minimal cou- 
pling between the matter lagrangian and the curvature 



((1 + z)h'^ - h^){l - ^hll^llllk) + (1 _ q^)^ + 



dh 
dz 

,n^{i + zfh 



'-{■ 



1 - go 
3 dh 



{l-qof ^l^zdz 
, 3(1 + zfn^f.h^ 



1 ,dh 2 d^h 
h dz dz"^ 

\3 



(1 



= a„(i + 2)ii + /l22) 



where subscript "a;" denotes differentiation with re- 
spect to variable x. The initial conditions are also taken 
in a way similar to the one described in the minimally 
coupled case in the previous section. The evolution of 




FIG. 5: Behaviour of the deceleration parameter q{z) as a 
function of redshift. The hues are for c = 2, 5, 10 from top to 
bottom, we assume ilm ~ 0.25 and go ~ —0.55. 



the Hubble parameter is shown in figure 3. Unhke the 
minimally coupled case, we do not have pathological be- 
haviour for any values of the parameter c. To compare 
the minimally coupled and non-minimally couple cases, 
we show the behaviour of the Hubble parameter for the 
two case for same values of c, fim and qq in figure 4. It 
shows that Hubble parameter in the non-minimally cou- 
pled case evolves slower than its counterpart in the min- 
imally coupled case. We also show the behaviour of the 
deceleration parameter for nonminimally coupled case in 
figure 5. Here also the universe is in an accelerating phase 
at present and smoothly joins the decelerating regime in 
the past. Here unlike the minimally coupled case, higher 
the value of the parameter c, the universe enters the ac- 
celeration regime earlier. 



together with the measurements of distance modulus by 
Type la supernova observations, to constrain our model. 
We use these two observational data to put constraints 
on our model parameter. In Figure 6, we show the la and 
2(7 contours in the qo — c plane with different choices of 
the density paremeter £7™. For minimally coupled case, 
there is always an upper bound for the parameter c as 
well as for the present day decceleration parameter qo- 
For example, with i7,„ = 0.25 the 2a upper bound for c 
is around 1.8. This upper bound shifts towards higher 
values as one increases flm- On the other hand, for non- 
minimally coupled case, things are different. Here there 
is an allowed range of c for every values of qo. For ex- 
ample, with rtm = 0.25 and qo — —0.9, the 2a allowed 
range for c is between 4.6 and 6.2. But as one increases 
f2m, this range shifts towards smaller value of c. Also for 
the minimally coupled case, our constraint on parameter 
c differs than that obtained by Ali et.al [2^ They solved 
the evolution equation in a different way than ours. They 
have assumed the universe is close to ACDM at higher 
redshifts, and put the initial conditions in the early time 
whereas we put the initial condition at present and solve 
it backwards. Also one of the initial conditions go is ac- 
tually one of the fitting parameters. Also they have used 
the slightly older Supernova data given by constitution 
set, which has 397 data points in comparison to our 557 
data points given by the Union2 set. Their analysis shows 
no constraint on c. at any level. It differs significantly 
from what we obtain. 



IV. CONCLUSION 



III. OBSERVATIONAL CONSTRAINT 

In this section we investigate the observational con- 
straints on our model parameters. We use the Supernova 
Type la data from the latest Union2 dataset consisting 
of 557 data point |23- The data consists of the distance 
modulus defined as fi = m — M = 5\ogdL + 25 where 
dL{a) is the luminosity distance defined as. 



dhia) 



dy 



y^H{yy 



(23) 



The other data we consider is the baryon acoustic os- 
cillations scale produced in the decoupling surface by 
the interplay between the pressure of the baryon-photon 
fiuid and gravity. For this we calculate the distance ratio 
Dy{z — 0.35)/£'„(z = 0.2) where Dy is given by 



Dy{z) 



dz' 



H{z) \Jo H{z') 



1/3 



(24) 



The SDSS observation gives Dy{z = Q.35)/Z)„(z = 
0.2) = 1.736 ± 0.065 [30i|. We use this measurements 



We have reinvestigated the f{R) gravity model where 
the curvature is minimally as well as in the case where 
gravity is non-minimally coupled with matter. We have 
assumed the Finder's Exponential form for /(i?) for our 
analysis. We have fixed the initial condition at present. 
We do not need any additional assumptions. By defini- 
tion, H/Ho = 1 at 2 = 0, where Ho is the present day 
Hubble parameter. The other initial condition fixes qo^ 
the present day deceleration parameter, and we take go 
as one of our fitting parameters. Hence our method of 
solving the evolution equation does not involve any ad- 
ditional assumption. First we check that for minimally 
coupled case, for c > 1.6, one can evolve the universe till 
any higher redshifts without any pathological behaviour. 
For lower values of c, there is some singular features in 
H{z) at different redshifts depending upon the parame- 
ter choices. For nonminimally coupled case, this changes 
and H{z) is regular till any higher redshifts for any choice 
of parameter values (we have taken f2„j between 0.25 and 
0.35, go between -0.9 to -0.55 and c between 1 to 50). For 
minimally coupled case, the constraint c > 1.6 to get reg- 
ular solutions for H{z) also satisfies the constraint com- 
ing from the local gravity tests [24j] . Also the evolution of 
the universe is as expected, showing accelerating universe 
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FIG. 6: la and 2cr contours in the qo — c plane for the minimally coupled (left panel) and non-minimally coupled (right panel) 
case using the observational data (explained in the text). The solid, dashed and dash-dotted lines are for flm ~ 0.25,0.3,0.35 
respectively. 



at late time, and decelerating universe in the past. We 
next use the the observational data coming from Type la 
supernova observations as well as the BAO peak measure- 
ments by SDSS. For supernova, we use the latest Union2 
compilation consisting 557 data points. For minimally 
coupled case, the constraints on c is completely different 
from what obtained by Ali et al. [2^ earlier. We have 
obtained a upper bound on c. This upper bound on c 
shifts towards the higher values as one increases ilm- For 
nonminimally coupled case, there is a range of allowed 
values for c, and this allowed range for c shifts towards 
for smaller values as one increases ilm ■ 



dia, through major research project grant (Grant No:33- 
28/2007(SR)). A.A.S also acknowledges the financial sup- 
port provided by the Abdus Salam International Center 
For Theoretical Physics, Trieste, Italy where part of the 
work has been done. ST and TRS acknowledge facili- 
ties provided by Inter University Center For Astronomy 
and Astrophysics, Pune, India through lUCAA Resource 
Center(IRC) at Department of Physics and Astronomy, 
University of Delhi, New Delhi, India. 



V. ACKNOWLEDGEMENT 

A.A.S acknoweledges the financial support privided 
by the University Grants Commission, Govt. Of In- 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 



S.Per lmutter et al, Astrophys. J. 517 (1999) 
[arXiv :astro-ph/9812133 .A. G. Riess et a i, 
tron.J.116 (1998) 1009 arXiv:astro-ph/9805201 



565 

As- 



i20T] 



D.N.Spergel et a/.[WMAP collaboratio n!, Astro - 

phys.J.Suppl. 170 (2007) 377 arXiv:astro-ph/0603449j 

U.Seljak et al [SDSS collboration],Phys.Rev. D 71 (2005) 

103515 arXiv:astro-ph/0407372 

D.J.Eisenstein et a/. [SDSS coUboration], Astro- [7] 

phys.J.633 (2005) 560 ar Xiv:astro -ph/0501171 

B.Jain and A.Taylor, Phys.Rev Lett. 91 (2003) 141302 [8] 

[arXiv:astro-ph/0306046 

E.J.Copeland, M.Sami and [9] 

S.Tsujikawa,Int.J.Mod.Phys.D 15, 1753 (2006): 



M.Sami, larXiv:0904.3445l V.Sahni and 

A.A.Starobinsky, Int.J.Mod.Phys.D 9, 373 (2000); 

T.Padmanabhan,Phys.Rep. 380, 235 (2003);E.V.Linder, 

asrto-ph/0704.2064; J.Frieman, M.Turner and Huterer, 

|arXiv:0803.0982 ; R.CaldweU and M.Kamionkowski, 

|arXiv:0903.0866; A.Silvestri and M.Trodden, 

arXiv:0904.0024 

C.Armendariz-Picon, T.Damour, and V.Mukhanov, 

Phys.Lett.B 458, 209 (1999) 

J.Garriga and V.F.Mukhanov, Phys.Lett.B 458, 219 

(1999) 

T.Chiba, T.Okabe, M.Yamaguchi, Phys.Rev.D 

62,023511 



[10] C.Armendariz-Picon, V.Mukhnov, and P.J.Steinhardt, 

Phys.Rev.Lett 85, 4438(2000) 
[11] C.Armendariz-Picon, V.Mukhnov, and P.J.Steinhardt, 

Phys.Rev.D 63,103510 (2001) 
[12] T.Chiba, Phys.Rev.D 66, 063514 (2002) 
[13] L.P.Chimento and A.Feinstein, Mod.Phys.Lett.A 19, 761 

(2004) 
[14] L.P.Chimento, Phys.Rev.D 69,123517 (2004) 
[15] R.J.Scherrer, Phys.Rev.Lett. 93, 011301 (2004) 
[16] A.Y.Kamenshchik, U.Moschella, and V.Pasquier, 

Phys.Lett.B 511, 265 (2001) 
[17] N.Bilic, G.B.Tupper, and R.D.ViolIier, Phys.Lett.B 535, 

17 (2002) 
[18] M.C.Bento, O.Bertolami, and A. A. Sen, Phys.Rev.D 66, 

043507 (2002) 
[19] A.Dev, J.S.Alcaniz, and D.Jain, Phys.Rev.D 67, 

023515(2003) 
[20] V.Gorini, A.Kamenshchik and U.Moschella, Phys.Rev.D 

67, 063509 (2003) 
[21] R.Bean and O.Dore, Phys.Rev.D 68,23515 (2003) 
[22] T.Muhamaki, M.Manera and E.Gaztanaga, Phys.Rev.D 

bf 69, 023004 (2004) 
[23] A.A,Sen and R.J.Scherrer, Phys.Rev.D 72, 063511 

(2005) 
[24] T.P.Sotiriou and V.Faraoni, arXiv:0805.1726 gr-qc]; 

S.Nojiri and S.D.Odinstov, arXiv:0807.0685 ; L. Amen- 

dola, R. Gannouji, D. Polarski and S. Tsujikawa, Phys. 

Rev. D 75 083504, (2007) ; W. Hu and L Sawicki, Phys. 

Rev. D 76 064004, (2007); A. A. Starobinsky, JETP Lett. 



[25] 



[26] 
[27] 

[28] 

[29] 
[30] 



86 157, (2007); S. Tsujikawa, Phys. Rev. D 77 023507, 
(2008) ; G. J. Olmo, Phys. Rev. D 72 083505,(2005); 
A. L. Erickcek, T. L. Smith and M. Kamionkowski, Phys. 
Rev. D 74 121501,(2006) ; V. Faraoni, Phys. Rev. D 74 
023529, (2006); T. Chiba, T. L. Smith and A. L. Erick- 
cek, Phys. Rev. D 75 124014,(2007); P. Brax, C. van de 
Bruck, A. C. Davis and D. J. Shaw, Phys. Rev. D 78 
104021,(2008); I. Thongkool, M. Sami, R. Gannouji and 
S. Jhingan, Phys. Rev. D 80 043523, (2009); R. Gan- 
nouji and D. Polarski, JCAP 0805 018, (2008); S No- 
jiri, S D Odintsov and D saez-Gomez, Phys. Lett. B 
681, 74,(2009); G Cognola, E Elizalde, S D Odintsov, 
P.Tretyakov and S Zerbini, Phys.Rev.D 79, 044001, 
(2009); S Nojiri and S D Odintsov, arXiv:0706.1378 
O.Bertolami and J. Paramos larXiv :0805.1241k2[gr-qc]: 
O.Bertolami, F. Lobo and J. Paramos, Phys. Rev. 
D, 78, 064036, (200 8); O. Bertolam i, J. Paramos, T. 
harko and F. Lobo, arXiv:081 1.2876, O. Bertoami and 
M.C.Sequeira, Phys. Rev.D, 79, 104010, (2009); V. 
Faraoni, Phys.Rev. D. 80, 124040 (2009). 
E.V.Linder Phys. Rev. D 80, 123528 (2009). 
Baojiu Li, and John. D Barrow, arXiv: 1007.3086 gr-qc]; 
J.D. Brown, Class. Quantum Grav. 10, 1579, (1993). 
Amna Ali, Radouane Gannouji, M.Sami and Anjan 
A. Sen, arXiv:1001.5384v2[astro-ph.CO]. 
R. Amanullah et a l., Astrophys.J. 7 16, 712-73 (2010). 
W.J.Percival et tt/, larXiv:0907.1660 astro-ph.CO] 



